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Abstract.  The  diocotron  spectrum  for  a simplified  model  of  Malmberg-Penning  traps  that  includes 
compression  effects  due  to  end  curvature  is  investigated  herein.  Performing  an  initial  value  treat- 
ment, we  find  that  there  is  a class  of  length  profiles  for  which  the  linearized  eigenvalue  equation 
of  the  model  can  be  integrated  in  quadratures  (integrable  case).  In  this  case,  there  is  only  algebraic 
growth  when  the  effective  angular  frequency  has  a maximum  (hollow  profile)  or  a minimum,  and 
the  model  is  mathematically  equivalent  to  the  zero  curvature  (2D  Euler)  case.  Furthermore,  we  study 
profiles  that  are  slightly  different  from  the  integrable  one  (the  difference  being  characterized  by  a 
small  parameter,  £),  finding  that  the  frequency  of  the  unstable  / = 1 mode  scales  as  e2/3.  Analytical 
calculations  and  numerical  simulations  are  found  in  remarkable  agreement. 


INTRODUCTION 

The  analogy  between  the  2D  Euler  equations  for  an  incompressible  and  inviscid  fluid 
and  the  classical  equations  for  the  evolution  of  a non-neutral  plasma  (2D  drift-Poisson 
model)  in  a Penning  trap  is  well  known  [1].  According  to  the  linear  theory,  the  2D  drift- 
Poisson  model  has  demonstrated  to  be  satisfactory  only  for  (exponentially  unstable) 
perturbations  with  azimuthal  mode  number  />  1.  On  the  other  hand,  it  can  be  proved 
theoretically  that  perturbations  with  mode  number  l = 1 are  always  stable,  regardless 
of  the  initial  equilibrium  density  profile  [1].  Experiments,  however,  show  that  the  linear 
growth  of  the  / = 1 mode  is  exponential  [2], 

New  fluid-dynamics  models  have  been  developed  for  a non-neutral  plasma,  trying  to 
solve  the  problem  of  the  / = 1 diocotron  instability  [3,  4].  These  works  show  that  a 
possible  explanation  of  the  instability  comes  from  the  finite  curvature  of  the  ends  of  the 
plasma  column  due  to  the  confining  voltage  (i.e.  compression  effects). 

The  present  study  investigates  the  modification  of  the  low-frequency  branch  of  the 
diocotron  spectrum  when  these  compression  and  finite  electron  temperature  effects  are 
taken  into  account. 
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PHYSICAL  MODEL 


We  focus  our  attention  on  a simplified  version  of  the  model  developed  by  some  of  the 
authors  in  Ref.  [4].  In  normalized  units,  the  model  is  the  following: 


<)a  „ 

Vi=e-xV1(J>eff 

->  O 

4>efr  = 0 + «log 


4)/= 0 


(1) 


For  a detailed  derivation  of  this  model  we  refer  to  Refs.  [4,  5,  6],  Compression  effects 
are  retained  in  the  terms  depending  on  the  normalized  temperature  a (assumed  to  be 
uniform)  and  on  the  normalized  effective  plasma  length  £o(r).  The  relation  between 
physical  and  dimensionless  quantities  can  be  found  in  Ref.  [5].  The  effective  plasma 
length  is  assumed  fixed  in  time  and  this  makes  the  model  2D. 

The  motivation  of  this  study  is  that  this  approximated  model  is  still  able  to  capture 
the  features  of  the  / = 1 diocotron  instability  (as  shown  in  Ref.  [4])  but,  since  it  is  2D, 
it  is  possible  to  perform  analytical  calculations  to  match  numerical  simulations.  Our 
goal  is  to  get  some  insights  on  the  effects  of  both  temperature  and  length  profile  on  the 
diocotron  spectrum  that  will  help  us  in  the  study  of  the  complete  model  of  Ref.  [4], 


ALGEBRAIC  INSTABILITY  IN  INTEGRABLE  CASE 


First,  we  perform  an  analysis  similar  to  the  one  by  Smith  and  Rosenbluth  [7].  Applying 
the  Laplace  transformation  to  the  linearized  time-evolution  equation  that  can  be  deduced 
from  model  (1),  one  obtains: 


(to  - IQ) 
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where  the  prime  means  derivative  with  respect  to  r and  8(f),,)  and  So  (r,0)  are  respectively 
the  Laplace  transformed  perturbation  of  the  potential  and  the  initial  perturbation  of  the 
line  integrated  charge  density.  A new  effective  angular  frequency  has  been  defined 


Q(r)=(o,(r)  + -^f\  (3) 

r Lo(i) 

where  o>£  ( r ) = (j)'^  ( r ) /r  is  the  E x B equilibrium  angular  frequency  in  the  central  plane 
of  the  trap.  In  the  limit  of  £o  =const,  Eq.  (2)  reduces  to  the  classical  Euler  case. 

After  some  algebra,  it  can  be  shown  that  under  a certain  condition  Eq.  (2)  can  be 
integrated  in  quadratures,  i.e.  8<)>(1,(r)  = r[to-£2(r)]  is  a solution  of  the  problem  for 
/ = 1.  We  refer  to  this  case  as  the  integrable  case.  The  condition  is  that  the  following 
equation  be  satisfied: 
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FIGURE  1.  Effective  length  profiles. 
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where  % = L^/Lq.  In  this  case,  provided  that  the  effective  angular  frequency  has  a 
stationary  point,  the  system  is  algebraically  unstable  and  the  mode  perturbation  grows 
asymptotically  proportionally  to  \ft  (as  for  the  classical  2D  Euler  case  [7]).  This  result 
has  been  confirmed  numerically. 

Given  the  equilibrium  initial  parameters,  Eq.  (4)  can  be  used  to  obtain  the  integrable 
length  profile  and  in  general  this  profile  will  be  different  from  the  one  that  comes  from 
the  equilibrium  solution  for  a given  configuration  of  the  Penning  trap.  In  order  to  make 
these  two  profiles  as  close  as  possible  inside  the  plasma  radius  (which  is  the  most 
important  region,  as  pointed  out  in  Ref.  [3])  we  use  the  boundary  conditions  x(0)  = 0 
and  x^O)  — Xocq(0)  /-G),eq(0)-  Note  however  that  one  could  in  principle  choose  the 
geometry  of  the  trap  (which  determines  the  equilibrium  profile)  to  match  the  integrable 
length  profile  XojntegW- 

Figure  1 shows  the  length  profiles.  The  equilibrum  length  profile  is  calculated  for 
the  same  geometrical  and  physical  parameters  of  Ref.  [4],  with  the  choice  of  confining 
voltage  V = 50  V and  temperature  0 = 7.5  V.  Accordingly,  we  choose  the  normalized 
temperature  a = 0.42  and  the  normalized  equilibrum  density  profile  [3,  4]: 
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with  no(0)  = 6.2,  fi  — 3 and  rp  = 0.59. 
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INITIAL  VALUE  TREATMENT 


In  Eq.  (2),  the  Laplace  transformation  of  the  perturbation  of  the  potential  is  defined 
as  S(J)0, (r)  = J'0+<*V<,)'S(|>(r,/)r//,  while  the  inverse  Laplace  transformation  is  given  by 

8<j>(r,r)  = where  r)  is  a large  enough  real  positive  number,  such 

that  all  irregular  points  of  the  function  8<i>(l)  are  located  below  the  line  Imo)  = q in  the 
complex  plane  of  to. 

We  introduce  the  function  t|/(r,o))  = r8<j>,„(r)  and  solve  the  second  order  equation  for 
induced  by  Eq.  (2).  Unlike  8<t>„,,  the  solutions  for  \{/  are  bounded  at  r - — y 0,  which  is 
necessary  for  Eq.  (6)  below. 

We  perform  our  analysis  in  the  vicinity  of  the  stationary  point  rs  (to*  = Q(rt)),  by 
using  the  method  of  matched  asymptotic  expansions  for  \j/(r,co).  For  simplicity,  the 
initial  perturbation  is  chosen  to  be  localized,  8o(r,0)  = a,8(r—  r,),  as  it  is  done  in  Ref. 
[8].  This  method  gives  the  following  result 


^t)  = -rrh)[ 
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. . , _ , 371 . / Q I, 

v|/r(r,(o)t|/s(/?H.,co)  - i 1 


3/2 


Vr  {R»;  0))\l/r(r,0>)  — V4/r(/?u  ,0))\[/v(r,  to)  X 


Vr(r„<rt)vi/,(0,G))  - \ps(r,,oi)vp,.(0,(D) 


Vr(0,CO)  f-T'(^yp-)  ' + - fi(0,(l))fr(^rifl) 


ck\) 

2n 


• (6) 


Here  the  contour  C goes  around  the  vertical  branch  cut  attached  to  the  point  to  = cov  in 
the  complex  plane  of  co.  \j/r (/-, co)  and  \|/s(r,(n)  arc  two  linear  independent  solutions  of 
the  homogeneous  part  of  the  Laplace  transformed  equation  such  that  t|/f  ~ 0{{r-rsf) 
and  t|tj  ~ 0(\/{r-  rs))  in  the  vicinity  of  rs.  We  have  also  introduced  A = to  - (i).s  and 
z/2£2| 


Q lr.= 


dr 2 


r=rs 


Two  merging  singular  points  are  present  when  r approaches  rs  and  the  contour  of 
integration  in  the  complex  plane  of  r should  always  go  between  these  two  points  [9]. 

We  expand  in  power  of  A the  integrand  of  Eq.  (6)  and  perform  the  integral.  Each  term 


of  the  expansion  contributes  as 


f Ave-'<0'  — - t 
Jr  2n 


-V-l  -l<Ojf 


. The  algebraic  instability 


(v  < —1)  is  only  possible  if  the  term  proportional  to  A~3/2  in  the  denominator  of 
expression  (6)  vanishes,  i.e.  if  t(/r(0,O))  = 0 or  \\ir (/?„,,  tn)  = 0.  In  the  integrable  case 
t(7r  = r(Q.(r)  — (Dj)  and,  therefore,  t|/r(0,o))  = 0.  Thus,  in  the  integrable  case,  the  leading 
term  of  the  expansion  has  V = -3/2  (from  the  numerator  of  Eq.  (6))  and  the  calculation 
of  the  contour  integral  leads  to  growth  « y/t.  Exact  calculations  of  the  contour  integral 
(6)  recover  the  asymptotic  formula  obtained  by  the  stationary  phase  method  in  Ref.  [7]. 

When  the  length  profile  deviates  only  a little  from  the  integrable  profile  (we  refer 
to  this  case  as  non-integrable  case),  £o  = (1  - e)A).imcg  + £A),cq>  where  e is  a small 
parameter,  the  perturbation  to  the  solutions  of  Eq.  (2)  are  « e.  In  particular,  8\jir  e and 
\|/r(0,O))  « £.  Since  t|/r(0,d))  ^ 0,  algebraic  growth  does  not  occur  in  the  non-integrable 
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Nr=2000 


FIGURE  2.  Frequency  shift  of  the  unstable  I = 1 mode,  log-log  scale. 


case  (continuum  modes  near  the  endpoint  to  = co^  phase  mix  according  to  t */2  in  this 
case).  The  zeroes  of  the  denominator  in  Eq.  (6)  determine  the  discrete  eigenfrequencies. 
One  can  see  that  for  small  e the  eigenfrequencies  exist  with  A « e2/3  (as  predicted  by 
Smith  in  Ref.  [10]  for  an  ad-hoc  model  and  seen  in  Ref.  [3]).  These  eigenfrequencies  are 
located  close  to  the  edge  of  the  continuum  cos. 

Numerical  simulations  have  been  done  to  compare  with  analytical  predictions.  Specif- 
ically, Fig.  2 shows  the  frequency  shift  of  the  real  frequency  (cor)  of  the  unstable  l = 1 
mode  with  respect  to  the  edge  of  the  continuum,  for  e between  0 and  1 . We  perform  the 
best  fit  with  the  least  squares  method  only  for  the  points  marked  with  triangles,  since  for 
very  small  e there  is  a loss  of  accuracy  and  for  large  E the  perturbation  theory  fails.  We 
find  C0j  — (Or « E0-70,  in  good  agreement  with  the  result  to*  — tor « e2/3  obtained  from  Eq. 
(6). 

Figure  3 shows  the  growth  rate  (y)  of  the  unstable  l — 1 mode  as  a function  of  E. 
Remarkably,  we  find  y^  E0  65,  in  excellent  agreement  with  y<*  e2/3.  One  can  also  notice 
that  the  integrable  case  is  a very  specific  case  and  as  soon  as  one  moves  a little  bit  from 
this  profile,  an  exponential  growth  is  found. 

Furthermore,  another  discrete  mode  appears  at  the  lower  edge  of  the  continuum  (in 
the  integrable  case  this  mode  is  the  stable  diocotron  mode).  As  predicted  analytically 
and  confirmed  numerically,  its  real  frequency  shift  (with  respect  to  the  lower  edge  of  the 
continuum)  scales  proportionally  to  £. 


REFERENCES 

1.  Davidson,  R.  C.,  An  Introduction  to  the  Plivsics  of  Nonneutral  Plasmas,  Addison-Wesley,  Redwood 
City,  1990. 

2.  Driscoll,  C.  F„  and  Fine,  K.  S.,  Phys.  Fluids  B 2,  1359-1366  (1990). 

3.  Finn.  J.  M„  del-Castillo-Negrete,  D.,  and  Barnes,  D.  C.,  Phys.  Plasmas  6,  3744-3758  (1999);  Finn, 
J.  M.,  del-Castillo-Negrete,  D.,  and  Barnes,  D.  C.,  Phys.  Rev.  Lett.  84.  2401-2404  (2000). 


333 


Nr=2000 


FIGURE  3.  Growth  rate  of  the  unstable  / = 1 mode,  log-log  scale. 
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